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Abstract: The nonperturbative effects in the high-energy processes involving strongly 
interacting particles are studied within the instanton liquid model of the QCD vacuum 
(ILM) by using the Wilson integral framework. The detailed analysis of nonperturbative 
contributions to the electromagnetic quark form factor is presented considering the struc- 
ture of the instanton induced effects in the evolution equation describing the high energy 



behaviour of the form factor. It is shown that the instantons yield in high energy limit 



the logarithmic corrections to the amplitudes which are exponentiated in small instanton 
density parameter. By using the Gaussian interpolation of the constrained instanton so- 
lution, we show that the all-order multi-instanton contribution is well approximated by 
the weak field limit result. The role of the instantons in high energy diffractive quark- 
quark scattering, in particular, in formation of the soft Pomeron, is also considered. We 
show that within the ILM the C— odd diffractive amplitude is suppressed as 1/s compared 
to the C— even one. The further applications of the developed approach in studying the 
nonperturbative effects in high energy hadronic processes are briefly discussed. 
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1. Introduction 

The powerful methods of the perturbative Quantum ChromoDynamics (pQCD) have been 
developed in order to describe the processes involving strongly interacting particles at high 
energies (for a review and comprehensive description of the methods, see, e.g., The 
total cross section of the e + e~ annihilation and the logarithmic violation of scaling in deep 
inelastic scattering became the classical tests of pQCD already in the lowest orders of 
expansion in the strong coupling constant and nowadays there are no doubts that the 
QCD Lagrangian provides a proper basis for a quantum field theory of strong interactions. 

At the same time, the present status of pQCD does not allow one to consider it as the 
only tool for investigation of the hadronic properties even at highest energies accessible at 
modern machines. The perturbative methods should be supplied by certain information 
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that can not be obtained directly from pQCD calculations. For instance, a nontrivial 
situation arises if, in order to make predictions of pQCD reliable, it is necessary to resum 
the soft part of the quark-gluon interaction to all orders. Moreover, in several situations 
the applicability of pQCD can be definitely justified only at asymptotically high energies, 
while in experimentally accessible region the nonperturbative effects are rather important 
and even dominant. 

Meanwhile, the intermediate energy regime — where Q 2 is larger than the typical hadronic 
scale determined by Aq£< d and lower than the characteristic scale of the chiral symmetry 
breaking: A.q CD ~ A-lonf < Q 2 < ^xSB' ^ s more convenient for detection of the non- 
perturbative phenomena. The complicated interplay of nonperturbative effects can lead 
in this regime to formation of the constituent quark which is, in a sense, an intermedi- 
ate object between color-neutral hadron and pointlike structureless partons, associated 
with the fundamental QCD particles — quarks and gluons. Indeed, the form factors of 
constituent quarks were recently extracted || from the JLab experiment data Q on the 
inelastic Nachtmann moments of the unpolarized proton structure function F^x^Q 2 ). 
And it was found that the size of a constituent quark determined from the data is about 
0.2 — 0.3 fm that corresponds to the mean instanton size in the instanton liquid model 
(ILM) H |?J which is also about 0.3 fm. In addition, the experimental and phenomenolog- 
ical investigations of the electromagnetic quark form factors at low and moderate energies 
can shed light on the problem of scaling violation in deep inelastic scattering processes. 

The requirement of certain nonperturbative supplement for explicit perturbative cal- 
culations even in high energy domain as well the study of the hadron processes at low 
and moderate energies appear to be a natural ground for development of nonperturbative 
methods. Coming down in the energy scale more and more powers of the strong coupling 
constant has to be taken into account. Moreover, in the intermediate energy region the 
power corrections come into play with coefficients that are very sensitive to the intrinsic 
hadron structure. Typically, the coefficients of the expansion in powers of the coupling 
constant and the inverse momentum transfer squared are the hadronic matrix elements of 
quark-gluon operators normalized at low energy scale and have to be found by nonpertur- 
bative methods. The dependence of these matrix elements on the energy scale is governed 
by the evolution equations that are determined within pQCD for different hard processes. 
These equations start to be applicable at momentum transfer squared of order 1 GeV 2 or 
higher where the strong coupling constant becomes small. So, it is necessary to find the 
initial data for the evolution equation that is essentially a nonperturbative problem. The 
important nonperturbative effects in this energy region are treated by different approaches: 
QCD sum rules, lattice QCD, quark models, etc. Therefore, the study of the role of the 
nonperturbative input in investigation of the processes with strongly interacting particles 
is not only an interesting theoretical problem, but is an important task for phenomenology 
of hadronic physics. 

It is naturally to relate the nonperturbative effects to the nontrivial structure of the 
QCD vacuum. In the last decades, a great progress has been made in study of the QCD 
ground state and a number of important results have been obtained that connect the 
properties of the vacuum with the hadron characteristics treating the QCD vacuum in the 
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framework of the instanton liquid model || |?J . Considering QCD vacuum as an ensemble of 
instantons, one can describe a number of the low-energy phenomena in strong interactions 
on qualitative and quantitative levels [||, || . The importance of the instanton induced effects 
in the strong interaction is also supported by lattice simulations ||, 10]. The instanton 
picture is generally considered as a fruitful and perspective framework for hadronic physics. 
The role of instantons in the hard hadronic processes has been studied intensively, both 
theoretically and experimentally. The contact with the perturbative QCD results becomes 
possible providing the unique information about the quark-gluon distribution functions 



in the QCD vacuum and hadrons at low energy normalization point 11, |12|, 13]. The 
perspectives for an unambiguous experimental detection of instanton effects are believed 
to be optimistic and promising. 

Although the QCD vacuum is known to play an important role in the high-energy 
scattering processes, the direct investigation of these effects remains a difficult task. The 
idea that the nontrivial vacuum structure could be relevant in high energy hadronic pro- 
cesses was explicitly formulated in the context of the soft Pomeron problem in [14, 15, 16], 



and further developed using the eikonal approximation and the Wilson integral formalism 



m 



In the present work, we report our recent results on investigation of the instanton 
induced effects in high energy regime for the electromagnetic (EM) quark form factor, 
i.e., the amplitude of the quark elastic scattering in an external color singlet gauge field. 
Then, we investigate the role of instantons in the diffractive quark-quark scattering and 
formation of the soft Pomeron. All the considered cases manifest a similar structure and 
are studied within the unified framework — the Wilson integral approach, which allows to 
study both perturbative and nonperturbative effects on the same ground. The method 
of path-ordered Wilson integrals is known as a powerful (and sometimes unique) tool in 
QCD which reformulates the theory in terms of the gauge invariant quantities — the Wilson 
loops — while the gauge fields are considered as chiral fields in the space of all possible loops 
|19| j. The Green functions, amplitudes, and cross sections can be expressed completely in 
terms of the Wilson integrals over contours with geometry determined by specific kinematics 



in an intrinsically non-diagrammatic, {i.e., nonperturbative) fashion |20j, 21, |22f|. 
2. Electromagnetic quark form factor 

The behavior of the hadronic form factors in various energy domains is one of the most 
important questions in the strong interaction phenomenology. The quark form factor has 
important phenomenological applications since it enters the cross sections of a number of 
high energy hadronic processes [||. For example, the total cross section of the Drell-Yan 
process (normalized to the deep inelastic one) is determined by the ratio of the time-like 
and space- like form factors [23, p3]: 



u n 
a PM 



F q(Q 2 ) 



2.1) 



F q (-Q< 

where <J^ M is the n-th moment of the cross section calculated within the parton model. 
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The electromagnetic quark form factors are determined via the elastic scattering am- 
plitude of a quark in an external EM field: 

M, = F q [(p! -p 2 ) 2 ] u(Pihf,v(P2) - G q [(px -pi?] u{p l ) <T ^ {Pl 2 ~ P2)u v{ P2 ) , (2.2) 

where u(p\), v{p2) are the spinors of outgoing and incoming quarks, and a^ u = [7^,7^/2. 
The kinematics of the process is described by two invariants (see Fig. la): 

s = (pi +P2? = 2m 2 (l + coshx) , t = -Q 2 = (pi - p 2 f = 2m 2 (l - cosh X ) , (2-3) 

where x ls the scattering angle. In this work we assume that both the momentum transfer 
—t and the total center-of-mass energy s are large compared to the quark mass: 

(P1P2) 2> Pi, 2 = m 2 , or coshx > 1 • (2.4) 

In this regime the Pauli form factor G q is power suppressed and will be neglected. However, 
it should be emphasized that in low and moderate energy domains it becomes important 
and there arise interesting perturbative and nonperturbative effects (see, e.g., recent works 

The EM quark form factor is one of the simplest and convenient objects for investi- 
gation of the double logarithmic behaviour of the amplitudes in QCD in the high energy 
regime. From the methodological point of view, it requires a perturbative resummation 
procedure beyond the standard renormalization group techniques. Besides this, the resum- 
mation methods developed for this particular case can be applied to many other processes 
which possess the logarithmic enhancements near the kinematic boundaries. Similar re- 
summation approach is also used in the study of the near-forward quark-quark scattering 



and the evaluation of the soft Pomeron properties [18|. In the latter case, the nonleading 
logarithmic terms become quite important. 

The first example of large logarithm resummation was given in QED by V.V. Sudakov 
for the case of an o/f-shell fermion in the external Abelian gauge field in the double- 
logarithmic approximation (DLA), where the terms of order of (a n ln 2n Q 2 ) are taken into 
account while the contributions from 0(a n ln 2n_1 Q 2 ) are neglected. The exponentiation 
of the leading double logarithmic result was found in |27f | . In DLA the exponentiated form 
factor behaves as a rapidly decreasing at high momentum transfer function. This means 
that the elastic scattering of a quark by a virtual photon is strongly suppressed. The 
exponentiation for the on-shell form factor in the Abelian case was obtained in the DLA 
in |^9| (the analytical calculation of the full electron vertex in two loops was performed in 

W)- 

These results eventually have been generalized for a more complicated non-Abelian 
gauge theory. First, the resumed DLA corrections in the QCD perturbative series were 
found to be consistent with the exponentiation in pH] (the inelastic on-shell form factor 



with emission of one and two gluons was calculated in the same context in [31]; the role of 
the quark Sudakov form factor in the description of e + e~ one-photon annihilation in quarks 
and gluons was considered in DLA in |p2| ), and the all-order non-Abelian exponentiation 
has been proved in [33|. 
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The question if the non-leading logarithmic terms could upset the DLA behaviour 
required a further work 34]. The all-(logarithmic)-order resummation was performed in 
the Abelian case and the exponentiation was demonstrated in [B^]. In |3(| the non-Abelian 
all-order exponentiation for the so-called hard part of the cm-shell form factor has been 
shown within the powerful factorization approach 1 . However, the status of the soft part, 
containing all the infrared (IR) and collinear singularities and, as a consequence, all possible 
nonperturbative effects, remained unclear. The important results on the IR properties of 
the QCD quark form factors was obtained in [37] within the Wilson loop approach. In 
these works, the soft part of the form factor has been presented as the vacuum averaged 
exponent of the path integral of a gauge field over the contour of a special form — an angle 
with sides of semi-infinite length (Fig. la). The use of the gauge and renormalization groups 
invariance allowed to derive the perturbative evolution equation describing the high energy 
behaviour of the form factor taking into account all (power unsuppressed) parts of the 
factorized amplitude, both for the on- |J8| and off-shell |3S!J cases. It was shown that the 
leading asymptotics is controlled by the cusp anomalous dimension which arises due to the 
multiplicative renormalization of the soft part, and can be calculated within the Wilson 
integral formalism up to the two-loop order |p7|] . Note that within the Wilson integral 
approach, the non-Abelian exponentiation can be proved independently [40], what is an 
important advantage of this framework. 

The efficiency of the Wilson integral approach has been successfully demonstrated in 



a series of works [22, f4l}| . In these papers the non-diagrammatic framework was developed 
that allows to calculate the fermionic Green's functions, Sudakov form factors, amplitudes 
and cross sections in QED and QCD completely in terms of the world-line integrals, and 
thus avoid complicated diagrammatic factorization analysis. 

The results presented above allow one to conclude that the leading high energy be- 
haviour of the quark form factor in a non-Abelian gauge theory is completely determined 
by the perturbative evolution equation, and is given by the fast decreasing exponent: 

2C F 



~ exp 

where 



0o 



lnQ z lnlnQ^ + 0{lnQ z ) 



(2.5) 



Nj-1 11N C - 2n f 
° F ~ ~2N7' 00 ~ 3 



This rapid fall off is not changed by any other logarithmic contributions [38], g9], 41]. 
However, the non-leading logarithmic corrections are nevertheless important for evaluation 
of the numerical value of the form factor at moderately large values of momentum transfer. 
Some of them are of a purely perturbative origin (e.g., sub-leading logarithmic terms), 
while the others can be attributed to the nonperturbative phenomena. 

In the present work, we try to advocate the point of view that the true nonpertur- 
bative effects can be taken into account consistently in the evolution equation, and there- 
fore they yield the non-vanishing subleading (perhaps, parameterically suppressed, but 

1 Note, that in this work the case where a time-like photon with large invariant mass decays into a 
quark-antiquark pair was considered, however it can be easily shown that the results remain true as well in 
the case of a quark scattering in an external EM field. 
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still logarithmic) contributions ~ In Q 2 to the high energy behaviour of the amplitudes. 
Further, we analyze another possible source of contributions which can be considered as 
"nonperturbative" — the IR renormalon ambiguities (there are plenty of papers on this sub- 
ject, for the most recent reviews see ^2|). We demonstrate explicitly that they produce 
the corrections with different IR structure compared to that one generated by instantons. 
Moreover, as we show below these renormalon effects disappear in the dimensional regu- 
larization |^4| and in the analytical perturbation theory p3[ , what means that they could 
be merely treated as artifacts of the incompleteness of the perturbative series resummation 
procedure. 

In the following Sections we describe the consequences of the RG invariance of the 
factorized form factor and derive the linear evolution equation considering the nonper- 
turbative input as the initial value for the perturbative evolution. We show that in the 
dilute phase the all-order single instanton contribution is exponentiated in small parameter 
of instanton density. Then, these nonperturbative effects are estimated in the weak-field 
approximation within the instanton model of QCD vacuum. By using the Gaussian sim- 
ulation of the instanton profile function we show that the weak field result approximates 
well the all-order multi-instanton contribution. The large-Q 2 behaviour of the form factor 
is analyzed taking into account the leading perturbative and instanton induced contribu- 
tions. The consequences of the IR renormalon ambiguities of the perturbative series and 
their relevance within the context of some analytization procedures are also studied. Fi- 
nally, the Wilson integral techniques is applied to evaluation of the instanton contributions 
to the high-energy behavior of the ^-scattering amplitude in the Regge regime, that is the 
soft Pomeron problem. 



3. Evolution equation for the quark form factor 



The classification of the diagrams with respect to the momenta carried by their internal 



lines allows one to express the form factor F q in the amplitude ( |2.2| ) in the factorized form 



m 36, m 



F q (q 2 )=F H (q 2 /fi 2 ,a s )-F s (q 2 /m 2 ^ 2 /X 2 ,a s )-Fj(fi 2 /X 2 ,a s ) , 



(3.1) 



where the hard, soft, and collinear (jet) part are separated by scale parameters /i 2 , A 2 . It 
is assumed the following relations between these scales q 2 » [i 2 » X 2 > 1. Note, that 
in the present paper all dimensional variables are assumed to be expressed in units of the 
QCD scale Aqcd, so that q 2 = Q 2 /Aq CD , etc. The arbitrary scale fj? is assumed to be 
equal to the UV normalization point. 

Within the eikonal approximation the resummation of all logarithmic terms coming 
from the soft gluon subprocesses allows us to express Fg in terms of the vacuum average 



of the gauge invariant path ordered Wilson integral |44|] 



F s (q 2 /m 2 , fJ 2 /X 2 ,a s ) = W(C X ; fj 2 / X 2 ,a s ) 



N r \ 



Vexp lig \ dx^A^x 



(3.2) 
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In Eq. ( |3,2[) the integration path corresponding to the considered process goes along the 
closed contour C x : the angle (cusp) with infinite sides (Fig. 1). We parameterize the 
integration path C x = {z^(t);t = [—00,00]} as follows 

*„(*) = | Vlt ' ~°° < * < °' (3.3) 
M ' \ vat, < t < 00. v ' 

The gauge field A^(x) = T a A%(x) (Tr[T a T b ] = \5 ab ) belongs to the Lie algebra of 
the gauge group SU(N C ), while the Wilson loop operator "Pexp (ig fdxA(x)) lies in its 
fundamental representation. 

The Wilson integral ( [3 .2D is multiplicatively renormalizable [[}?], 20]: 

W{C x -^ 2 /X\a s (i?)) = Z cusp (C x ;fL 2 /fi 2 ,a s (^ 2 )) ■ W bare (C x ; ft 2 / X 2 ,a s (fi 2 )) , (3.4) 

where p? is the UV cutoff, /i 2 is the normalization point, and A 2 is the IR cutoff. The 
presence of the IR divergence in ( |3.4| ) is a common feature of on-shell amplitudes in massless 
QCD. Therefore, we can define the cusp anomalous dimension T cusp : 

i T cusp (C x ;a s (fx 2 )) = -fi 2 -^ In W(C X ; fi 2 /\ 2 , a s {ix 2 )) = (3.5) 
= -^ 2 -^2^Z cusp (C x ,fl 2 /fi 2 ,a s (fi 2 )) . 



It can be shown that the cusp anomalous dimension (|3.5| ) is linear in the scattering angle 
X to all orders of perturbation theory in the large-g 2 regime [37|: 

T cusp (C x ;a s ) = Inq 2 T cusp (a s ) + 0(ln% 2 ) , (3.6) 



00 

'a. 



where 

r cnsp (a s ) = ^2 (~) CnAi , (3.7) 
7r 

C n are maximally non-Abelian color factors (C n = CfN™' 1 in lowest orders), and A n are 
some numerical factors. The observation ( p.6[ ) is crucial for the efficiency of the analysis 
of high energy behavior within the Wilson integral formalism, since calculating the cusp 
anomalous dimension only in leading (one-loop) order one obtains the dominant contribu- 
tion to the high energy asymptotics. 

The total form factor F q is renormalization invariant quantity and satisfies to the 
renorm group (RG) equation: 

/, 2 A F g (^,a s (» 2 )) = L*^+p(a a )^L) Fq ^ 2 , a^ 2 )) = , (3.8) 

which in the large-g 2 regime leads to the following relations 

2 d \ d]nF H ] _ 2 d \ dlnF s ] _ 1 . . . . 

J < 2 ^lno 2 ~ ^ du 2 dhxa 2 ~ 2 Lcusp[as) - ^ 



dfi 2 
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To avoid the problems with additional light-cone singularities in the soft part, we work 
with logarithmic derivatives in q 2 45]. The collinear part Fj being independent on q 2 
does not contribute to these equations. 



Then, from the Eqs. (3.S, 3.5 , p^q ) one finds after simple calculations |p8| : 
din F H (q 2 ) 



3 In q 2 

dlnF s (q 2 ) 
din q 2 



/A 4 (jC 
2 ^T cusp (a s (Z))+T(a s (q 2 )) , 



cusp 



>.(0) + 



dlnW np (q 2 
din q 2 



where the "integration constant" of the hard part reads 

T(a s (g 2 )) 



2 din F H (q 2 ^ 



dlnq 2 



and W np arises as the initial value of the soft part: 



(3.10) 
(3.11) 

(3.12) 



dlnW np (q 2 ) _ dlnF s (q 2 



dlnq 2 



din q 2 



(3.13) 



H 2 =\ 2 

and is the only quantity where, according to our suggestion, the nonperturbative effects 



take place [46, 47]. 



Then we get the q -dependence of the total form factor at large q : 



In 



F g (q 2 ) 



ln q -T. 



cusp 



q 2 f q o d£ 
In— / — T cusp {a s (t])) +ln 



W np (g 2 ) 
W n Jqfr 



(3.14) 

In the next subsection, we explicitly calculate the perturbative quantities entering Eq. 



(3.14) in the one loop perturbative QCD approximation and in the other parts we concen- 
trate on the nonperturbative part within ILM. 



4. Perturbative contributions to the Wilson integral 

The analysis of the hard contributions [36, ^] at large q 2 yields: 
d In F H (q 2 /fi 2 ,a s ) 



d In q 2 



-?H ln H)+ o( "* 



(4.1) 



For the hard "integration constant" ( |3.12| ) one has: 



rK(^)) 



2 ... _ 3a s (g 2 ) 



4 7T 



C F . 



(4.2) 
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The expression (|4.l| ) is IR-safe, while the low-energy information is accumulated in the soft 
part of the quark form factor Fg. The Wilson integral ( |3.2| ) can be presented as a series: 



W(C X ) =i + ^ (0 |Efe)7 I ~f 

c n _2 ^x x JO x 



where the function 6(x) orders the color matrices along the integration contour. 
The leading order of the expansion (|4.3| ) is given by expression (see Fig. lb-c): 



(4.3) 



9 



9_ 

' 2 A",. 
2 1 



1 -TT(T a T b ) f dxJ dy v V%{x-y) 



dx, 



dy v T>fj, v (x - y) 



(4.4) 



2 2N, 

where the gauge field propagator V^ u (z) in n-dimensional space-time (n = 4 — 2e) reads: 



ZV(*) 



T4(z)^(0) 



(4.5) 



WW e (C x ; e, /I 2 /A 2 ) = exp K^ e (C x ; e, /i 2 /A 2 ) + O(o*) 



The exponentiation theorem for non-Abelian path-ordered Wilson integrals [40] allows us 
to express (to one-loop accuracy) the Wilson integral (|3.2| ) as the exponentiated one-loop 
term of the series ( |4.3| ): 



(4.6) 



In general, the expression (|4.4j ) contains UV and IR divergences, that can be multiplica- 
tively renormalized in a consistent way [2C]. In the present work, we use the dimensional 
regularization for the UV singularities, and define the "gluon mass" A 2 as the IR regulator. 

The dimensionally regularized free propagator in covariant gauge reads (n = 4 — 
2e) , e > 0: 



4— n 



d n k 
(2vr)" 



-ikz 



*v,w - (i - oa 2 + »o] 



l fc 2 - A 2 + iO k 2 - A 2 + iO 

where £ is a gauge fixing parameter. It is convenient to use the representation 

D^(z) =g^d 2 A 1 (e,z 2 ,fi 2 /\ 2 )-d^d u A 2 (e,z 2 ,[l 2 /\ 2 ) , 



(4.7) 



(4.8) 



where p? is a parameter of dimensional regularization. 

By using integrals ( 13. ip - ( |13.6 ) from Appendix A, the dimensionally regularized 
formula for the leading order (LO) term ( |4.4P can be written as pcf : 

(4.9) 



w£l(C x ;e,fi 2 /\ 2 ,a s ) = 8na s C F h( X ){l - e)A 1 (e, 0, /x 2 /A 2 ) 
where h(x) is the universal cusp factor: 

Hx) = Xcothx - 1 , 



(4.10) 
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which at large-g 2 is given by: 

-.2 



lim h(x) — > X ^ l n ~~n ■ (4-11) 

X— >oo 777/ 



In Eq. fl4.9Q, for the perturbative gauge field one has 



A 1 ( £ ,0,AVA 2 ) = - T ^f44V (4-12) 



16vr 2 V A 2 / 1 - e ' 

The independence of the expression ( |4.9| ) of the function A2 is a direct consequence of the 
gauge invar iance. 

Thus, in the one-loop approximation one gets 



W bare (C x ; e, fl 2 /X 2 , a.) = l- ^C F h( X ) ( - - 7 £ + ln4^ + In ^ ) , (4.13) 



and the cusp dependent renormalization constant, within the MS'-scheme which fixes the 
UV normalization point, reads: 

Z cusp {C x -e,ii 2 /^\a s (ii 2 )) = l + ^^C F h{x)(^ . (4.14) 



Using the Eq. ( f4.9[ ), one finds the known one- loop result for the perturbative field, 
which contains the dependence on the UV normalization point /x 2 and IR cutoff A 2 (see 
e.g., mM): 



2 



W^\c x ^ 2 /\\a s {^)) = -'^^C F h{x)\n^ + 0{a 2 s ) . (4.15) 



Therefore, in the leading order the kinematic dependence of the expression ( ^4.4[) is factor- 
ized into the cusp factor h(x) |H - 

From the one-loop result ( |4.15 ), the cusp anomalous dimension which satisfies the RG 
equation (|3.5| ) in one- loop order is given by 2 : 

r^ sp M//)) = ^^c F . (4.16) 

7T 

Substituting into the Eq. Q3.14 ) the anomalous dimension ( |4. 16| ) with the strong 
coupling constant given in the one-loop approximation, one finds 



2C F 




[ Po 





exp L -_ ^ - l) - + ln, 2 (l -ln^|)] + W np{q 2 )] 

(4.17) 

Note, that the exponent in Eq. ( [4.17 ) has an unphysical singularity at A 2 = 1 (in dimen- 
sional notations, A 2 = Aq CD ), i. e., where the one-loop coupling constant a s (A 2 ) has the 

2 The cusp anomalous dimension is known up to two-loops, see |37[. 
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Landau pole. This feature can be treated in terms of IR renormalon ambiguities (see next 
Section), and is often considered as a signal of nonperturbative physics. In the present 
paper, we will consistently separate the sources of nonperturbative effects which can be 
attributed to uncertainties of resummation of the perturbative series from the "true" non- 
perturbative phenomena. An important example of the latter is provided by instanton 
induced effects within the instanton model of QCD vacuum, which is considered in the 
Section 6. 



5. IR renormalon induced effects 



As it was pointed out at the end of the previous Section, the perturbative evolution equation 
( 4.17 ) possesses an unphysical singularity at the point A 2 = 1. Therefore, it is instructive 



to study the consequences of this feature. It is known that the presence of the Landau 
pole in the one-loop expression for the coupling constant leads to the IR renormalons p2[ 
resulting in power suppressed corrections. In the present situation one can expect the 
corrections proportional to the powers of both scales: /i 2 and A 2 . We will treat here the 
power // 2 -terms to be strongly suppressed in large-g 2 regime, and focus on the power A 2 - 
corrections. To find them, let us consider the perturbative function Aj(e, 0, /2 2 /A 2 ) in the 
Eq. (|4.9|). The insertion of the fermion bubble 1-chain to the one-loop order expression 
(|4.4| ) is equivalent to replacement of the frozen coupling constant g 2 by the running one 



g —* g (k ) = 4ira s (k ) [48] (for convenience, we work here in Euclidean space. 



A.M^/A 2 ) = -4^/ ^^) f { ^ 2) . (5.1) 

By using the integral representation for the one-loop running coupling a s (k 2 ) = * dcr(l/k 2 ) ah , 
b = /?o/4-7r, we find: 

*, 2 ^2n 1 f<P?Yf°°, r(l - x - e)r(l + x + e) ( l\ x . , 

Ai(e,0,/2 2 /A 2 = -— : / dx — — ; r -hro • (5-2) 

H '' M/ ; A)(1-£)V / Jo {x + e)T(l-e) \\ 2 J v ' 

To define properly the integral in r. h. s. of Eq.(fT^), one needs to specify a prescription to 
go around the poles, which are at the points x n = n, n € N. Thus, the result of integration 
depends on this prescription giving an ambiguity proportional to (l/A 2 ) n for each pole. 
Then, the IR renormalons produce the power corrections to the one-loop perturbative 



result, which we assume to exponentiate with the latter [48|. Extracting from (5.2) the 
UV singular part in vicinity of the origin x = 0, we divide the integration interval [0, oo] in 
two parts: [0,5] and [6, oo], where S < 1. This procedure allows us to evaluate separately 
the ultraviolet and the renormalon- induced pieces. For the ultraviolet piece, we apply the 
expansion of the integrand in Ai in powers of small x and replace the ratio of T-functions 
by exp(-7£e): 



„ 1 ( hi 47T — JE + hi T2 ) r s / 1 \ x 



k,n=0 
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which after subtraction of the poles in the MS'-scheme becomes: 



Ar(0,/x 2 /A 2 ) 



1 



A)(l 



£ 

n=l 



In 



A 2 



n-l 



(5.4) 



In analogy with results of [49], this expression may be rewritten in a closed form as 



(fx 



-a; In 



-zlnA 2 



Then, using the relation 



one finds 
.2 d 



cw(% 2 ) 

9 In g 2 



2C>(1 



^r(o,/i 2 /A 2 ) 



9(7 y Sing 2 



exp 



4tt 



/3 a s (/i 2 ) 



(5.5) 



(5.6) 



(5.7) 



The exponent in the last equation yields the power suppressed terms (l/q 2 ) S in large-g 2 
regime. In LLA the Eq. (|5.6|) is reduced to: 



dWW {q 2 
dlnq 2 



2C F / In// 2 
In 



In A 2 



(5.8) 



The last expression obviously satisfies the perturbative evolution equation ( 4.17| ) . 

The remaining integral in Eq. ( |5.2| ) over the interval [5, oo] is evaluated at e = 
since there are no UV singularities. The resulting expression does not depend on the 
normalization point /i, and thus it is determined by the IR region including nonperturbative 
effects. It contains the renormalon ambiguities due to different prescriptions in going 
around the poles x n in the Borel plane which yields the power corrections to the quark 
form factor. 

After the substitution fi 2 = q 2 and integration, we find in LLA (for comparison, see 
Eq. (pTD): 



F r ™(q 2 ) = exp 



2C F 



Inq 2 (lnhV-1) - lnq 2 <S> ren (\ 2 ) 



(5.9) 



where the function $ ren (A ) = Ylk=a 0fc(VA ) accumulates the effects of the IR renor- 
malons. The coefficients cpk cannot be calculated in perturbation theory and are often 
treated as "the minimal set" of nonperturbative parameters. It is worth noting that the 
logarithmic g 2 -dependence of the renormalon corrections in the large-q 2 regime is factor- 
ized, and thus the Eq. ( |5.9| ) corresponds to the structure of nonperturbative contributions 
found in the one- loop evolution equation (|4 .17 ). On the other hand, it can be shown that 
careful account of partially resumed perturbative series yields, sometimes, cancellation of 
the leading power corrections associated with the leading renormalon contributions [48]. 
As the corresponding nonperturbative terms are calculated independently (e.g., by means 
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of the ILM) their direct relation to the IR renormalon ambiguities should be questioned. 
In our point of view, it allows us to separate the true nonperturbative (e.g., instanton 
induced, but not only) effects from that ones related to ambiguities of the resumed pertur- 
bative series. 

The latter conclusion can be illustrated by considering the consequences of an analyti- 
zation (AN) of the strong coupling constant |S| in the perturbative evolution equation. In 
this approach, the one- loop strong coupling constant a s ([i 2 ) is replaced by the expression 
which is analytical at fx 2 = 1 (i.e., at A.q CD in dimensional variables): 



The direct substitution of ( |5.10| ) into the evolution equation ( 3.14j ) yields (for brevity, we 
assume q 2 , = A 2 ): 



( Ini^V) = In q 2 ( In _ 1 )- * ln >E* + (5.11) 



In q 2 + I) (ln l) X l~^ ) ~ \ (1^ ? ~ In 2 A 2 ) - Li 2 (l - q 2 ) + Li 2 (l - A 2 ) 



The functions Li 2 in the resulting expression accumulate the power corrections of q 2 and 
IR scale A 2 , but does not exhibit a singularity at A 2 = 1. Therefore, it gives no room for 
IR renormalons ambiguities, at least in the considered approximation. Nevertheless, the 
power corrections of nonperturbative origin do contribute to the large-g 2 behaviour. Note, 
that the consequences of the analytization of the strong coupling constant in the IR region 
have been studied also in [5(| . 

Another possible way to avoid the Landau pole at the integration path was developed 
within the dimensional regularization (DR) scheme ]^4j]. In this case, the running coupling 
reads 

afW) = MOT ; -i] ' <5 ' 12) 

and for complex e , Re e < 0, it has the Landau pole at the complex value of /i 2 , thus this 
singularity moves out of the integration contour. In the limit e — > 0, the form factor reads 



1 24] (for comparison, see Eq. ( |4.17|) : 

if V) = 

exp - — - I + ln e fin q 2 - - J + ln q 2 (ln ln q 2 - l) - - ln ln q 2 J +0(e,e ln e) 

(5.13) 

This expression also leaves no room for any renormalon effects. At the same time, the 



stant" W np , ( 3.11 ), which is not directly related to the analytical properties of the coupling 
constant. 
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6. Instanton induced contribution to the Wilson integral 



In this Section we consider explicitly the nonperturbative effects within the instanton model 
of QCD vacuum. The instanton effects in the high energy QCD processes has been actively 
studied since the seventies |5l], ^] . Recently, the investigation of these effects was renewed 
with promising perspectives |26|, [h], 47, 53|, 54, ^5], 56, 57, 58, 59, |(|. The Wilson integral 



formalism is considered as a useful and convenient tool in the instanton applications mainly 
due to significant simplification in the path integral calculations if an explicit form of the 
gauge field is known. Another important feature of this approach is the possibility to make 
a correct analytical continuation of the results obtained in the Euclidean space (where the 
instantons are only determined) to the physical Minkowski space-time where the scattering 
processes actually take place. Namely, one maps the scattering angle, %, to the Euclidean 



space angle, 7, by analytical continuation |6T] 



A^x;p)=A°(x;p) — = -R^VjUx - *oWO - z ; p), (6.2) 



X^ij, (6.4) 

and performs the inverse transformation to the Minkowski space-time in the final expres- 
sions in order to restore the g 2 -dependence. 

Let us consider the instanton contribution to the function W np (q 2 ) from Eq.( 3.13j ). 
The instanton field has general form 

2 9 

where <fi(x) is the gauge dependent profile function, M a ' ) is the color orientation matrix, 
a a, s are the Pauli matrices, 77 „„ = e^ afll/ =F {l/2)e a b c ei, cfiu are 't Hooft symbols, and (±) 
corresponds to the instanton or antiinstanton solution 3 . 

The averaging of the Wilson operator over the nonperturbative vacuum is performed 
by the integration over the coordinate of the instanton center zq, the color orientation and 
the instanton size p. The measure for the averaging over the instanton ensemble reads 
dl = dM d 4 ZQ dn(p), where dR refers to the averaging over color orientation, and dn(p) 
depends on the choice of the instanton size distribution. Taking into account ( |6,2| ), we 
write the Wilson integral ([T^) with contour ([T^) in the single instanton approximation in 
the form: 

wi{x) = ■^Tr(0\exp[ia a (n a 1 a{vi,z ) + n%a{v 2 ,z ))}\0) = (6.3) 
= -^Tr(0|cos(a(fi,z ))cos(a(u 2 ,zo)) (I x I) - 
- g (n\ ■ n 2 ) sin(a(vi, z ))sin(a(v2, z )) (a a x a a ) |0), 



where (i = 1,2) 



\ 1) ±a U v 

7 sV nu V i 

S{Vi,Zo) 



^lu<4 , (6-4) 



3 Below we always consider only topologically neutral instanton (I) and antiinstanton (A) configurations, 
but since there are no differences in contributions of I and A in the processes considered both solutions will 
be called as instanton. 
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poo 

a(vi,z ) = s(vi,z ) dXip [(vi\ + (-If z ) 2 ; p c ] , 
Jo 



and 



S 2 (Vi,Z Q ) = Zq - (ViZoY 



(6.5) 



(6.6) 



where (vyu^) = cosh^ in Minkowski geometry. We omit the path ordering operator V in 
( |6.3| ) because the instanton field (|6.2|) is a hedgehog in color space, and so it locks the color 
orientation by space coordinates. Let us note that due to nonperturbative factor g^ 1 in 
the instanton field Q6.2| ) the Wilson integral ( |6.3[ ) is independent on the coupling constant. 

At this stage we need to average over all possible ways of embedding SU C (2) into 
SU C (3) by using the correspondence relations |32| 



((I X 1) SU M))SU(3) ~ q 



(I X I) 



sup) ' 32 



'5(7(3) 



(6.7) 



3 

{{a a x (J a )su{2))su{3) = g (A A x A A ) S{/(3) 



where A" 4 are generators of SU(3) algebra (A = 1,2, ...,8). Then we get 
1/4 

wj(x) = 3 Tr \° gCos(a(«i,«d))cos(a(u 2 ,2o)) (I x I) + 



1 

+ - 



-cos(at(«i, 2b))cos(a(u 2 , 20)) - sin(a(ui, ^))sin(at(v 2 , ^n"^ 



(6.8 

(A A x A A ) |o) 



The resulting gauge invariant contribution to the Wilson loop of the single instanton taken 
in all orders in gauge field becomes [46| 



wi{x) = 3 / dn (p) [ w l(x) + wi{x) - w J c (0) - w{ 



w cix) = J d 4 ZQ cos a(vi,zo)cos 0(^2,^0) > 

w lix) = ~ d 4 z (n1h%)sin a(v 1 ,z Q )sin a(v 2 ,z ) , 



(6.9) 
(6.10) 



where the normalized color correlation factor is 



VilA^>2^0 

s(vi,z )s(v 2 ,z )' 



(6.11) 



The Eq. ( |6.9| ) takes into account the subtraction of the infinite self-energy parts of the 
quark form factor. Thus, to calculate the instanton contribution one need to consider only 
the vertex corrections. 



7. Exponentiation of the instanton contributions in the dilute regime 



On the basis of the exponentiation theorem [ 40 1 for the non-Abelian path-ordered expo- 
nentials it is well known that perturbative corrections to the Sudakov form factor are 
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exponentiated to high orders in the QCD coupling constant. The theorem states that the 
contour average Wp{C) can be expressed as 



W P {C) = exp 



E(v)" E Cn(W)F n (W) 



n=l 11 W€W(n) 



(7.1) 



where summation in the exponential is over all diagrams W of the set W(n) of the two- 
particle irreducible contour averages of nth order of the perturbative expansion. The 
coefficients C n (W) oc CpN™~ 1 are the "maximally non-Abelian" parts of the color factor 
corresponding to the contribution coming from a diagram W to the total expression ( |7,1| ) 
in the contour gauge 4 , and the factor F n (W) is the contour integral presented in the 
expression for W. This means that the essential diagrams are only those, which do not 
contain the lower-order contributions as subgraphs and, as a result, the higher-order terms 
are non-Abelian. 

Let us now demonstrate how the single instanton contribution is exponentiated in the 
small instanton density parameter, treating the instanton vacuum as a dilute medium p4[ . 
The gauge field is taken to be the sum of individual instanton fields in the singular gauge, 
with their centers at the points zfs. In this gauge, the instanton fields fall off rapidly at 
infinity, so the instantons may be considered individually in their effect on the Wilson loop. 
Moreover, the contribution of infinitely distant parts of the contour may be neglected and 
only those instantons will influence the loop integral, which occupy regions of space-time 
intersecting with the quark trajectories. Since the parameterization of the loop integral 
along rays of the angle plays the role of the proper time, a time-ordered series of instantons 
arises and has an effect on the Wilson loop. Thus, the contribution of n instantons to the 
loop integral Wj(x) can be written in the dilute approximation as 

Wi n) ( x ) = Tr (u 1 U 2 ...U n U n l..U 2 ^U 1 ^ , (7.2) 
where the ordered line integrals Ui's 

U 3 {x) = T jexp ^ig J da v^A^a - Zj) \ exp (ig J dr v^A^(v 2 t - Zj)\ X 

are associated with individual instantons with the positions Zj's. Because of the wide 
separation of the instantons in the dilute phase and rapid fall off of fields in the singular 
gauge, the upper and lower limits of the line integrals are extended to infinity. The line 
integrals U l ^s take into account the infinitely distant part of the contour that goes from 
+oo back to — oo and in the singular gauge they are IP' = 1. For the integral is 

taken over the increasing (decreasing) time piece of the loop. 

Then, the expression is further simplified by averaging over the gauge orientations 
of instantons. The averaging is reduced to substitution of W by gjU^gJ 1 , where gj is 



4 The contour gauge is defined as a gauge where the condition Pexp [ig dz ll A /1 (z)) — 1 is fulfilled 
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an element of colour group, and independent integration of each gj over the properly 
normalized group measure is performed. Under this averaging one gets 

Unjjn^ ^ w (n) (%) = ^ ^jnjjn^ , ( 7 . 3 ) 

which is just the single instanton contribution as it is given by Eqs. (|6.9|). If the averaging 
is done in the inverse order, from n down to 1, the entire loop integral becomes the product 
of traces 

n 

W\ n \ x )^ lim n»f'(x). (7.4) 

n— >oo 

3=1 

Since the individual instantons are considered to be decoupled in the dilute medium, the 
total multiple instanton contribution to the vacuum average of the Wilson operator simply 
exponentiates the all-order single instanton term wi{\) in ( |6.9[ ), and one has 

Wi(x) = lim \l + -Mx)\ = exp[w/(x)]. (7.5) 

n— »oo in) 

Thus, we prove that in the dilute regime the full instanton contribution to the quark 
form factor is given by the exponent of the all-order single instanton result (see Fig. le). 
The exponentiation arises due to taking into account the multi-instanton configurations. 
As it is well known, in QED there occurs the exponentiation of the one-loop result due to 
Abelian character of the theory. In the instanton case, the analogous result arises since 
instantons belong to the SU (2) subgroup of the SU (3) color group and the path-ordered 
exponents coincide with the ordinary ones. 

The following comments are in order. First, the nonperturbative exponentiated ex- 
pressions are strictly correct only as long as the instanton density n(p) is small. Second, it 
is supposed that U(zq) is evaluated using the singular gauge form of A™ st . On the other 
hand, Tr (UU'j is identically the ordered loop integral for a single instanton and is gauge 
invariant. It is therefore legal to use the nonsingular gauge of A l ™ st in evaluating the trace 
(a more handing gauge for computation). 

8. Large-g 2 behaviour of the instanton contribution in the weak field ap- 
proximation 

Let us first consider the weak-field approximation to the instanton contribution to the 
quark form factor. In this limit the leading instanton induced term in Eq. ( |6.9| ) is given 
by general expression ( [4.4D (Fig. |2jc) with the corresponding dimensionally regularized 
instanton correlator 

= A ™" 4 / d <P) J (0r p)K(-k; P )e-* kz . (8.1) 

By using the Fourier transform of the instanton field embedded into the SU{3) color group 

K(k-P) = --R aa v ±a ^'(k 2 ;p) , (8.2) 
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where <p(k 2 ;p) is the Fourier transform of the instanton profile function (p(z 2 ;p) and 
(p'(k 2 ; p) is its derivative with respect to k 2 , and the properties of the color rotation matrices 
and 't Hooft symbols 

rt-r = 5^ , 7^< A = g^gpx - g^xgvp + e^ P x , 
one can express the instanton correlator (|8.1| ) in the form similar to Eq. (fh8]) 

V U Z ) = ^2 (9^d 2 ~ W A d z2 ) , (8-3) 

where 

Adz 2 ) = -A™ -4 J dn(p) [ 2 £(k 2 ; p)] 2 . (8.4) 

Then, applying the same method as described in perturbative case, one gets the LO in- 
stanton contribution in the form 



™r\x) = Wc h(x) J dn(p) Aj(0,p 2 X 2 ) , (8.5) 

where we use the same IR cutoff A 2 , while the UV divergences do not appear at all due to 
the finite instanton size. In the singular gauge with the profile function given by 

P 2 

WW = 2 1 2 l 2V ( 8 - 6 ) 

ar (x z + p z ) 

one gets: 

A 1 ^p 2 X 2 ) = ^[Hp 2 X 2 )^{p 2 \ 2 ) + ^{p 2 X 2 )] , (8.7) 



where 



1 f 1 dz 



/A 4 J z(l - z) 

and 



„2 \2\ 



A 2 ^0 



- ~p 2 X 2 (xz + y{l-z))T Mz+v{1 _ 



$ x (p 2 A 2 ) = V / dxdydz ~ r - ' — i^eP-A-iM-w-*;] i im ^(p^) = 

^— ; Jo n!n a 2 ->o 

n=l u 

(8.9) 

are the IR-finite expressions. 

At high energy the instanton contribution is reduced to the form: 



dlnW^q 2 ) vr 2 



J dn{p) p 4 [ln(p 2 X 2 ) $ (p 2 A 2 ) + $i(/9 2 A 2 )] = -^(A 2 ) . (8.10) 



ding 2 AN C 

Here we used the exponentiation of the single-instanton result in a dilute instanton ensemble 
(see [46] and the previous Section) and took only the LO term of the weak-field expansion 
( |4.4j ): = Wj + (higher order terms). 

In order to estimate the magnitude of the instanton effect we take the instanton size 
distribution found by 't Hooft [65] and multiply it by the exponential suppressing factor 
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(which was suggested in [ 67 ] in the framework of constrained instanton model and assumed 
in |68t| in order to describe the lattice data lit 



dn (p) = -r c n c 
P b 



2tt 



a s (fj, r 



2N C 



cxp 



2tt 



a s (p r 



(p/i r ) /3 6Xp (-27TC7P 2 ) , 



.11 



where the constant C Nc = 4.6/vr 2 exp(-1.679JV c )/ [{N c - 1)!(JV C - 2)!] « 0.0015, a « 
(0.44 GeV^) 2 is the string tension, (3 = 0q + 0(a s (fj, r )), and is the normalization point 
|6£ 



Given the distribution ( |3.11 ) the main parameters of the instanton liquid model, the 
instanton density n c and the mean instanton size p c , will read: 



/ dn(p) 
Jo 



C N T(f3/2 - 2) 



2tt 



2N C 



A, 



QCD 



V2 



ira 



(2vra) 2 , 



/ °°p dn(p) _ T((3/2 - 3/2) 1 



J °°dn(p) rCS/2-2) V2 



7TC 



1.12) 



U3) 



In Eq. ( |8 . 1 3| ) we choose, for convenience, the normalization scale /i r of order of the instanton 
inverse mean size p~ l , taking into account that the distribution function ( p. 1 1 ) is the RG- 
invariant quantity up to 0{ot 2 s ) terms 6S]. Note, that these quantities correspond to the 
mean size p c and density n c of instantons used in the model [p|, where the size distribution 



{.11 ) is approximated by the delta- function: 

dn(p) = n c 5(p - p c )dp . 



Thus, we find the leading instanton contribution ( 8.1C ) in the form: 

A 2 



where 



Bi(X 



K 



> _ K f ln 2 ^ 



2N r 



A 2 



l + O 



r(A,/2)[T(A)/2-2)f' 



2na 



0.74 , 



2 [r(A)/2 - 3/2)]4 
and we used the one loop expression for the running coupling constant 

2tt 



(3 ln p c A 

In ( 8. 15| ) we introduces the packing fraction parameter 

t 2 4 



(8.14) 

(8.15) 
(8.16) 

(8.17) 
(8.18) 



that characterizes diluteness of the instanton liquid. Within the conventional picture its 
value is estimated to be 

/w0.12, (8.19) 



if one takes the model parameters as (see §): 

n c 1/m" 4 , p c 1/3/m. 



.20) 
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Therefore, from Eqs. ( |8.10| ) and (5.9), we find the expression for the quark form 
factor at large-g 2 with the one-loop perturbative contribution and the nonperturbative 
contribution (the function W np in Eq. Q4.17 )) found in the instanton model: 



F q (q 2 ) = exp 



^C-F , 2i i 2 

hi q in In q 

A) 



Inq 2 (Bj 



2C F 



+ 0(\u lng 2 



F (q 2 ;X 2 



.21) 



Numerically the coefficient Bj is small factor as compared to the perturbative term 



Bi « 0.03 <C 



2C F 



0.24, 



(8.22) 



where to get this estimate we take the IR cutoff parameter A of order of Aqcd ~ 350 
MeV. It is clear, that while the asymptotic "double-logarithmic" behaviour is controlled 
by the perturbative cusp anomalous dimension, the leading instanton correction results 
in a finite renormalization of the subleading perturbative term (Figj3|). Note, that the 
instanton correction has the opposite sign compared to the perturbative logarithmic term. 

It is also important to note that the results obtained are quite sensitive to the way one 
makes the integration over instanton sizes finite. For example, if one used the sharp cutoff 
then the instanton would produce strongly suppressed power corrections like cx (l/q)^°. 
However, we think that the distribution function (8.11) should be considered as more 



realistic, since it reflects more properly the structure of the instanton ensemble modeling 
the QCD vacuum. Indeed, this shape of distribution was recently advocated in [|^, 67] 
and supported by the lattice calculations 



(for comparison, see, however, [7C, 71]). 



9. All-order calculations of the Wilson loop for Gaussian profile 

The weak field limit used in the previous Section may deviate from the exact result. Here, 
we are going to test its accuracy considering the properties of the single instanton contri- 
bution to the Wilson loop, w/(x)> defined in Eq. ( |6.9| ), which contains contributions from 
the gauge field taken in all orders. 

We have to note that in realistic instanton vacuum model there are two essential effects: 
stabilization of the instanton density with respect to unbounded expansion of instantons in 
size (see Eqs. (|8TD) and <KW>), and screening of instantons by surrounding background 



fields. To take into account these features we approximate first the narrow instanton size 
distribution by the 5- function as in Eq. (|1|) and assume that the integration over the 
instanton size is already performed. As it was discussed above, the screening effect modifies 
the instanton shape at large distances leading to the constrained instantons |67|. To take 
into account this screening and to have also simpler analytical form for u>/(x), we use i n 
this Section the Gaussian Anzatz for the instanton profile function 



Vg{x 2 ) = \e- x2 /£ . (9.1) 



P? 



The parameters in this expression are fixed from the requirement to reproduce the vacuum 
average value of the (^(0)yl"(0)) which is finite for the instanton field in the singular 
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gauge and equals 5 : 

(^(0)^(0)) = I2^pln c . (9.2) 

Below we equal the parameter p c to unity and restore the dependence on it at the end of 
calculations. 

The phases fl6,5| ) corresponding to the profile ([O]) are calculated as 

a (v h2 , zo) = s(v lt2 , Zo )^e- s2 ^^eifc [ T (v h2 z )} , (9.3) 
where with definitions ( |3.3| ) one has 6 

(viz ) = 24, (v 2 z ) = z 4 cosh x + IZ3 sinh x, (9.4) 
s 2 (v 1 ,z ) = z\ + z\, s 2 (v 2 ,z ) = (z 3 coshx - iz 4 sinhx) 2 + z\ . (9.5) 

We take the reference frame where the scattering point is in the origin and define the 
scattering vectors of the quark as 

vx = (1, 0, 0_l) , v 2 = (cosh x,i sinh x,0_l) , (9.6) 

v \ = v 2 = 1 ! (vi«2) = coshx , 

where the velocities t>i = pi/m and t>2 = p 2 /m determine the scattering plane. In this 
system the instanton position is given by zq = (z^, z%, z±). 

However, due to exponentially large oscillations at large x that occur during integra- 
tion over instanton position, it is not easy to use the closed form for the phases (|6.5| ) in 
general case and for the Gaussian profile ( |9.1| ), in particular. We have to note that these 
complications do not arise in the case of high-energy elastic quark-quark scattering consid- 
ered in |)3| . In the latter case, the angle dependence simply factorizes from the integrand 
in the asymptotic regime (see the next Section). 

In order to cancel exponentially large oscillations, we need to integrate first over the 
instanton position zq. For this purpose, by using the explicit form (|9.1|) we expand the 



expressions (6.10) in powers of the phases ( |6.5| ) and change the order of integrations: 
w cKX) 2^ L > (2n)!(2ro)! ' s{X) ^ ^ K ' (2n + 1)! (2m + 1)! ' 

n=lm=l \ / \ I n=0m=Q V ) \ ) 



(9.7) 



where 



^2n^2m' 
Oil a 2 



(*) = J II***] II *K J d^(v 1 ,z ) S ^(v 2 ,z )e-i {v ^- Z0) + ^ +z ^ \ 

(9.8) 



i=i -« j=i 



5 This average is infinite in the regular gauges. Thus, the singular gauge realizes the minimal Landau 
gauge for the instantons. 

6 In these expressions and below we make formally analytical continuation to the Minkowski space that 
is inverse to the transformation (3.1). 
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and 



2n+l 2m+l\ 



a 



(x) 



/ IT ^ / TT dA i / ^O^" ("1' ^) S 2m (v 2 , Z ) si 2 (z ) 

Jo H Jo fj, 3 J 



(«lA,-2:o) 2 + (^2A^+2:o) 2 ] 



(9.9) 



Note, that the integrands possess the symmetry with respect to the change of variables 

z 3 ~^ z 3 cosri X + i%4 sinh x , Z\ — > z 4 cosh x — IZ3 sinh % • (9.10) 

Above we introduced the notation for the color spin correlation factor (its normalized 
definition given above in (6.11)) 



Si 2 (zo) = Z3 (2:3 coshx — 1Z4, sinhx) + z\ coshx • 



(9.11) 



Let us further make the change of variables A by introducing the total and partial 
lengths 



N 



JV 



i=l 



with new measures given by 



dLL 



N-l 



L 



d{x} 



{K}m 



M 



V 



(9.12) 



N 1 



dL'L' M - 1 



d{y}u , 



(9.13) 



where J d{x}w = Jq dx\... Xl "' XN 2 dx^-i- 

After these transformations we come to the general expression for the phase averages 
(|9~^ ) and d]| with arbitrary M and A 



/■OO /-00 

(ofaf)(x)=/ <&/ dL'e-TOr U ^^Giv^lI.i'.xj^Mli)^^^) , (9-14) 

JO JO 



with definitions 

Fn,m{L) 



(9.15) 



1 J d{x}x exp < 



1 



A + M 



A' 



v 



L 2 + (A + M - 2) >H ~ 4$3 A * A 3 



G Nj m{L, L', X ) = J d 4 z s 2n ( Vl ,z ) s 2m (v 2 , z ) sll(z )e-( N+M ^ +z ? +z ?) , (9.16) 

where Z4 = 24 — L ~^_^^ h x and z 3 = 23 + ^jyqr^p ■ The definition for Fm,n(L') is similar 
to one of Fjv,m(-^)) the power 77 in Gn,m(L, L' , x) is equal to one for (x) and to zero for 
wi (x)- 

In principle, the integral in Gn,m(L, L' ,x) may be done analytically for the Gaussian 
profile in arbitrary order of expansion for any finite values of x- I* 1 this case, after the 
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d zq integration there remain no complex valued variables anymore and the Euclidean 
(transverse) and Minkowski (longitudinal in x) dependencies get factorized. But in practice 
only a few first terms can be analyzed. Nevertheless, as we will see below, in the large x 
limit it is possible to make the partial summation of the double sums in the expressions 
Q9.7| ). Let us demonstrate these statements in more detail. 

Consider the lowest order one-loop and two-loop contributions given by (Fig. ||) 

w\x) = ^ (-(aia 2 )( X )|s + g<a?a 2 >(x)|s + ^(«l«l)(x)|s + \{a$a%)(x)\s + •• 

(9.17) 

where (af 1 )(x)|s = ( a x a 2)(x) ~~ ( a i a 2 The simplest one- loop diagram (Fig. 

^|a) corresponds to the first term in ( 9.17 ) ( N = M = 1 term in Eq. ( |9.14 )). In that case, 



the integrals may be done analytically (by using the Eqs. ( |HD1) and (HJ)) and the final 
result reduces to the weak field expression (|8.5| ) with the function Ai given by 



4 

Thus, in the lowest order we have the result (/ = 7r 2 n c p^,) 



Af (x 2 ) = -fA e - 2 l2 P l (918) 



g fl) (x) = -^Hx) , (9.19) 



12 



and its small and large x behavior 



4 1,1} (X - 0) = -^x 2 , (X - oo) = -L x . (0.20) 



It is important to note, that due to screening of instantons at large distances the result 



(9.18) is IR-finite contrary to the case of the pure instanton solution in (| 

If one neglected the color spin factors in (9.16), one would get another dependence on 
the scattering angle 

4^ (X) = ~{ ( 1 - , (9-21) 



6 \ sinh^ 
(1.1) t n\ f 2 (1,1) / \ / 



thus getting the asymptotic which is not enhanced by the large factor x- 

The rest of Eq. ( p,17| ) corresponds to the two-loop diagrams shown in Fig. 
(N = 1, Af = 3 in ( |9l4D ) and in Fig. |c (N = M = 2 in flOD). Then, the functions 
Fn,m{L) and Gn,m{L,L' ,x) becomes 

W = ex P (~\ L2 ) - = yfexp (-I//)erf(A 

Jo •/ 
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G 2 , 2 (L,L', X ) = J d A z s 2 (v 1 ,z )s 2 (v 2 ,z )e-<^ +t ' 2+z '^ = (9.23) 
vr 2 

= — j^" \L' 2 sinh 2 x{L 2 sinh 2 x — 6) — 8LL' sinh 2 x cosh x — 6£ 2 sinh 2 % + 52 + 8 cosh 2 x] , 
Gi,3(L,^,x) = / ^o S 2 (^i^o)^ 2 (^o)e- 4 W + '' 2+ ^) = 

7T 2 

= — [-LL' 3 sinh 4 x + 10L' sinh 2 x(£' cosh x + L) - 60 cosh x] • 

We compare the full result and the one-loop and two-loop approximations to it at small x 
in Fig. ||. 

In Fig. H we also present the results corresponding to the calculations without color 
spin factors, s (vi,zo). In the latter case, the coordinate integral J d^zo in ( |9.16| ) may be 
performed easily with the result 

G' NtM (L,L', X )= {M + N) 2 ■ ( 9 - 24 ) 

The next order calculations may be done in a similar way. Finally, it may be shown that in 
the limit of large scattering angle x the asymptotics reads w Gs (x — * oo) ~ const. Thus, in 
this case one has weaker asymptotics than the asymptotics with color spin factor included. 

In the following, we are interested in the limit x ~* 00 ) where the coefficient of x 
is free of the light-cone singularities and therefore it has a well defined limit for the on- 
shell quark momenta, p 2 = p| = 0- To find this asymptotics, we take into account that 
Fn,m(L — > 0) ^ and Gn,m(L, L') are polynomials in L and L'. With these properties 
and in the large x limit, the L and L' integrations may be performed analytically 

POO POO 

lim / dL dL 1 (LL') n e~ aLL ' cosh *F N M {L)F M n(L') = 
x-*°°Jo Jo 

= 7 jr^F NtM (p)F M jf(P) , ( 9 - 25 ) 

[a cosn x) 

roc roo 

lim / dL dL'(LL') n L m e- aLL ' coshx F NM (L)FMN(L') = 
x^°°Jo Jo 

= 7 '—^ttFm,n(0) / dLL m Fn,m{L) . (9.26) 

(a cosh x) J 

The important result is that the only diagonal terms with equal powers of L and V provide 
the leading in x asymptotics. 

Let us consider the contribution to the asymptotics from the terms with m = and 
arbitrary n of w{. (x) (Fig. 0a). This contribution is reduced to the element 

POO /'OO P 

(al n+1 a l 2 ) = -coshx / dLL 2n / dL 1 / d{x} 2n+1 e~ L2 ^ X "~ L ' 2 • (9.27) 
Jo Jo J 

■J d A Z {zl-iZzZ A t a nhx + zl) ( z | + 4) n e -[2(n+lK 2 -2^(L-L'co S h X ) + 2 i23 L' S inh X ] 

The zq integration gives (see Eqs. ( |13.7| ) and ( |13.8|) ) 

o , , „ f°° „ f°° f , ,„ L 2 + L ,2 -2LL' coshx 

(a 2 ^ +1 a\) = -7T 2 coshx / dLL / dL' / d{x} 2n+1 e- h2 ^ ^~ L \ . 

Jo Jo J 
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fe! 



k + 1 / iL' sinh x 



+ 



+ 



k =o [2(n + l)l" 

^'sinh X N 2 ^- fc+1 ) n ^- +1 



2(n + l) V 2(n+l) 



2(n-Jfc) n-fe 



2(n + l) 



zLtanhx /iL'sinhx 



-^2(ra-fc+l) 



s=0 

2(n-fc)+l n-Jfe 



s=0 

iL' sinh x 
2(n + l) 



-^2{n-k) 



2(n-fc+l) 



s=0 



V 



2(n + l) V 2 ( re + !) 
where we introduce notation 



2(n + 1) 
(il/ sinhx) z 



(9.28) 



(9.29) 



The Eq. ( |9.28| ) is still a general expression. To find its large- x limit, we need to analyze, 
according to (|9.25| ), the coefficients of the maximal powers of the diagonal terms (LL') 2n . 
In the first and fourth terms in the brackets of ( |9.28| ), the only terms with k = s = 
provide the leading x asymptotics and the other one gives the subleading contribution. 
The second and third terms in the brackets give dominant asymptotics if k + s < 1. Their 
sum provides the leading x-asymptotics, while the terms of higher powers in xcoshx are 
canceled. Keeping only leading terms in the Eq. (9.28), one gets 

7r 2 cosh x 



I 2n+l 1\ 



X^oo 



16(n + 



^3 J™ dLL 2n J™ dL' J d{x} 2n+l e- L2 ^ 



2(n + l) 



-X 



2n 



2n + 3 



LV sinh 2 (x) 
(n + 1) cosh(x) 



2(n + l) 

As x ~~ * 00 j all integrals can be done analytically by using Eq. ( |9.25| ): 

(_l)n+l 7r 2( 2n )! 



-X 



(9.30) 



(9.31) 



ix^oo 2 4ri+3 (n + l) 2 

Now, let us consider the contribution to the asymptotics of w^, (x) from the terms with 
m = 1 and arbitrary n (Fig. |7|b): 

poo pOO f f'l 

* 2n a 2 ) = cosh 2 x / dLL 2n - 1 / dL'L' / d{x} 2n / dye' 1 * E^^'V+M 2 ) . 
Jo Jo J Jo 



«1 



d 4 zo 



1 + zf + 2iz^Z4 tanh x~ z 4 tanh 2 x 



cosh x 



7t 

• ^ (^n z 2A: z 2(n--fc) e _[2(n+l)zg+224(i-L' cosh x)+2iz 3 sinh x] 32) 
fe=0 

By using the table integrals Eq.( [13.7 ) and ( 13. 8j ), the Gaussian integration over zq yields 

pOO pOO /• f'l 

{a\ n a 2 2 ) = vr 2 / dL dL' {LL') L 2 ^ / d{x} 2n / dy e ^ 2 . 



•e 



£ ^ +£ /z_ 2LL f coshx 
2(n+l) 



fc! 



fc=0 



[2(n + l)] 



ra+3 



L u sinrr x 
2(n + l) 



(n-k) 
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-^2{n-k) s ( L sinh x 3 



-fc-i 

C^" fc (n — k — s)y s+1 LL' sinh 2 x cosh % 



(9.33) 



s=0 

with y given by Eq.( |9,29| ). Again, the large x limit of this expression stems from the 
diagonal terms (LL') 2n 

/ 2n 2\ I 

L 2 + L 12 -ILL' cosh x 

. e 2(^+T) [n(n + l)(2n - 1) - 4n(n + 1)LL' sinh x + (LL' sinh x) ] • (9.34) 

After integration one finds that the coefficient of the leading asymptotics is equal to zero, 
and therefore 

(aj n al)\ ^ = const . (9.35) 

Moreover, it is possible to show that further leading asymptotic terms appear only if n > 3 
and m > 3, but they are highly suppressed numerically. 



Thus from fl6.9|) , (|9.7|) and (|9.3l[) , we find that the leading correction to the quark form 
factor is given by 



oo 



WG(q 2 ) = -Y ^ cn/ | x - jrlng 2 = -1.0053 lng 2 , (9.36) 



^ 16 n (n + l) 2 (2n + l) 12 ^ '12 

n=0 v ' K ' 



that means that the weak field limit ( |9.20| ) is a good approximation for the Gaussian profile 
function ( |9.1[ ). This estimate of the logarithmic coefficient — Bq ~ 0.01, is compatible with 
the estimate ( |8.22| ) obtained in the single instanton approximation. 

We have to comment that in the general (non-Gaussian) case the weak-field limit may 
deviate from the exact result. Nevertheless, it is reasonable to expect that by taking the 
instanton solution in the singular gauge which concentrate the field at small distances and 
allows us to prove the exponentiation theorem for the Wilson loop in the instanton field 
|46], one gets a good numerical estimate of the full effect. Thus, the resulting diminishing of 
the instanton contributions with respect to the perturbative result appears to be reasonable 
output. The analysis of all-order instanton contribution performed in the last part of this 
Section for a Gaussian profile function shows that the weak field approximation can be 
justified, but an additional investigation of this problem is required. 



10. Instanton model of Pomeron 

Soft hadronic collisions are described successfully within the Regge phenomenology, with 
the Pomeron exchange being dominating at high energy. The Pomeron is considered as 
an effective exchange in the t channel by the object with vacuum quantum numbers and 
with positive charge parity C = +1. That is why the idea that the nontrivial structure 
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of the QCD vacuum is relevant in describing its mechanism. To illustrate this idea, let us 
consider the high energy diffractive quark-quark scattering, where there is a hope that for 
small momentum transfer the nonperturbative effects give dominant contribution. One of 
the simplest models of the Pomeron is based on the use of exchange by two nonperturbative 
gluons. The nonperturbative part of the gluon propagator is given by (in the Landau gauge) 



: A a Jx)A b v (0) : 



1 5 ab 



g 2 N*-lJ (2tt) 



d A k 



-ikx 



k^k 1 



G 



2 ) • (io.i) 



In the Abelian gauge model considered originally by Landshoff and Nachtmann ||16|| , 
the nonperturbative gluon propagator (without a color factor) G np {k 2 ) is related to the 
correlation function describing the gauge invariant gluon field strength correlator (nonlocal 
gluon condensate). In general non- Abelian case this correlator has the form 



: Gfj, u (x)V exp 



ig I dz a A a (z) 
o 



G pa (0) : 



/d 4 k 
-—^- lkx [{D (k 2 ) + D^k 2 )) k 2 (g w g vry - g^g up ) + 
(2vr) 

-\-D\(k ) [kpkpg vu — kpk a g V p + k v k a gp p — k^kpg^^)^ , 



(10.2) 



where the first tensor structure is called non- Abelian part and the second one is Abelian 
part. Indeed, in the Abelian gauge model without monopoles Do(k 2 ) = 0, and D\{k 2 ) = 



G np (k ). It is this property that has been used in []T6| to relate the Pomeron properties to 
the value of the gluon condensate. 

However, in the non- Abelian model one has the opposite situation. Really, for the 
QCD instantons we find [66, 67] Di(k 2 ) = and Do(k 2 ) is nonzero. In the realistic 
model of the QCD vacuum, where the interaction with vacuum fields of large scale, R, 
is important, the instanton ceases to be exact solution of the equations of motion, but 
the so-called constrained instanton approximate solution (CI) can be constructed It 
has been shown that the constrained instanton has exponentially decreasing asymptotics 
at large distances (~ R). The constrained instanton has topological number ±1 like an 
instanton; however, it is not self-dual field. Thus, in the realistic QCD, the small non-zero 
part of D\{k 2 ) appears. Very similar results have been found in the lattice simulations of 



the gluon field strength correlator |72 ]. 

Thus, within the non- Abelian models there is no direct connection of the gluon prop- 
agator to the gluon field strength correlator. So, let us explicitly consider the instanton 
part of the gluon propagator. The Fourier transform of the instanton field is defined as 



where 



~ Air 2 i f 00 

( j ) (k 2 ) = — r dxx 3 J 2 (\k\x)ip(x 2 ) 



K Jo 



(10.3) 



(10.4) 
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and J2{z) is the Bessel function. The function 4>(k 2 ) is related to the function ip(k 2 ) defined 
in dp) by 



-2i- 



dk 2 



(10.5) 



The explicit form of the Fourier transform of the pure instanton solution is well known (in 
the singular gauge) 



)=i 



(4tt) 
fc 4 



{Pckf 



K 2 (p c k) 



4>\k 2 ) 



t(47r) 2 
fc 4 ' 



A: 2 



0. 



oo . 



The (constrained) instanton profile may be chosen in the form (c.f. Eq. (|B.G|)) 



Phi 



xr 



x 2 (x 2 + p 2 JI (x 2 )) 



, Phi* 2 ) = 2 ( 



(10.6) 



(10.7) 



where K%(z) is the modified Bessel function. The constrained solution saves its form at 
short distances, but changes it at large ones: 



0ci(k 2 ) = 

where the constant Ici is given by 

Ici 



l 4-R A I C i , k 2 ^0, 



fc 4 ' 



k 2 — > oo 



du u 2 ipci{uR 2 ) . 



(10.? 



(10.9) 



Now the Fourier transform of the single instanton contribution to the gluon propagator 



( |10.1[) becomes 

G np (k 2 ) = -An c k 2 ^ 2 {k 2 ) . (10.10) 

Thus, we see that the nonperturbative gluon propagator ( 10.1| ), (10.10) and gluon field 
strength correlator ( 10.2 ) are quite different functions, and the relation between them that 
is valid in the Abelian gauge model is destroyed in the non-Abelian case. 

From ( |10.6 ) and ( 10. 8| ) one gets the asymptotics of the instanton part of the gluon 
propagator 



G 



I ( r 2\ 



{{2v) A n cP ik- 2 , 
\ (4^) 4 n c fe" 6 , 



k 2 ^0 



■ G CI (k 2 ) 

k — > oo | ( i~) ' n, .k " . /.r — > oo 

(10.11) 

Calculating (in a very similar way as in the Landshoff-Nachtmann model) the invariant 
T— matrix element of the quark-quark scattering at large energy, s, and small transferred 
momentum, t, we get 



7r 4 n c R 8 j2 £.2 



16 ± CI I% > 

(4vr) 4 n c A:- 6 , 







k 2 



(<lk(P3)qi(P4) \T\ q m (Pl)qn(P2)) 
with 



il(t) u(p 3 )j ll u(pi)u(pA)j fJ -u(p2)SkmSi n , (10.12) 



dk± 



G 



2^ 



(10.13) 
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where 

s = (jp x +p 2 ) 2 = 2m 2 (l + coshx), t = -(p 3 - pi) 2 = -q 2 ± , 

G(k 2 ) is defined in ( |10. 10 ) with k 2 — > k\. Except of the numerical coefficient, this ex- 
pression is in agreement with the Landshoff-Nachtmann formula. This agreement is due 
to the specific features of the instanton induced interaction. From the infrared behavior 
of the instanton induced propagator (10.11), it is clear that 1(0) is infinite for the pure 
instanton solution ( 10.6j ), but it is finite for the constrained instanton solution ( 10.7] ) . This 
fact, also noted recently in [68], was one of the arguments to construct the constrained 
instanton solution. The form of I(t) is presented in Fig. ^ for the constrained instanton 
configuration. 

It is possible to show that the result ( |10.12| ) is the weak field limit of the more general 
expression. Indeed, the quark-quark scattering amplitude may be expressed in terms of 
the vacuum average of the gauge invariant path ordered Wilson integral [17, [13] 



T% n (s,t) = -2is f d 2 b ± e^W^ n ( 7 ,b 2 ± ) , 



where the Wilson line function WM n is given by 



(7,63.) 







Vexp lig I dx^A^x) 



Co 



ki 



(10.14) 



(10.15) 



In Eq. ( 10.15| ), the corresponding integration path goes along the closed contour C qq : two 
infinite lines separated by the transverse distance b± and having relative scattering angle 
X- We parameterize the integration path C x = {z^(A); A = [—00, 00]} as follows 



v 2 \ + b± 



-00 < A < 00 , 
— 00 < A < 00 



(10.16) 



with (t>ii>2) = coshx and b± being the impact parameter. 

By making the steps similar to the previous Section, one arrives to the expression for 
the Wilson line function (c./. Eq. ( |6.10D ) (see Fig. ||) 

"1 



W( 7 ,bl) =n c {^ c ( 7 ,6i)(lxl) + i 



-w c {l,b ± ) + w s (-y,b A 



w s ( r 7,b 2 L ) 



d zq [cos a(v\, zq) — 1] [cos a(v2, Zq) — 1] , 



d zo(n"n2)sin a(vi,zo)sin 0(^2,-^0) 



(X A x X A ) } (10.17) 
(10.18) 



where the color correlation factor is 

(viv 2 )(z ,z - b±) - (v 1 z )(v 2 z ) 



h\h a 2 



s(vi,zq)s(v 2 ,z ) 



The phases are defined as 

/oo 
dX (p [(zq + viX) 2 ;p] , 
-00 

/•oo 

a(v 2 , z ) = s(v 2 , z ) / dX (p [(z - v 2 X - b±) 2 ; p] 



(10.19) 

(10.20) 
(10.21) 
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with 

s 2 (vi,z ) = zl - {vxZq); s 2 (v 2 ,z ) = (z - 6_l) 2 - (v 2 z Q ). 

By means of the change of variables (214COS7 — Z3sin7) — > 24, the energy dependence is 
trivially factorized 



1 



sni7 



-io c (7 = 7r/2, b 2 ± ) , w a (rf, b ±) cot 7 Ws{l = 7T/2, 



(10.22) 



At 7 = 7r/2 the above definitions are reduced to 



s 2 0i,3 ) 
n?n^( 7 = 
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tt/2) 


= -1 


+ z 2 , s 2 (t;2,2;o)- 
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a(ui,2d) 
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= Si 


I dX <p [z 2 + z\ 

i-00 


+ X 2 ;p] , 




-> 02 


= S2 


/ dX tp [zl + (z 

J —00 


-b) 2 ± + X 2 ;p] 



(10.23) 

The differential cross section of the quark-quark scattering is expressed through the 

d<Jn 



amplitude ( 10.17 ) as 



1 1 



dt 9 s 2 



•qq 



££|^(m)I 2 . 



(10.24) 

kl mn 

By inserting ( 10.17| ) into ( 10.24| ) and making analytical continuation to Minkowski space 
one finds 



da qq {t) 2 2 



coth 2 x F 2 (t) + *^X Fe (t)F.<t) + T-T^W 
3 smh x o sinh x 



, (10.25) 



where 



F a (t) = j d 2 b ± e ib ^w s ( 7 = tt/2,6 2 ) , F c (t) = j d 2 b ± e ib ^ 



w c (j = 7r/2,b 2 ± ) . (10.26) 

In the asymptotic limit (sinhx ~ s , cothx — > 1) the result ( |10.25| ) coincides with that one 
found in [53 

§«§-WM. (10-27) 

In the weak field limit we reproduce the one-loop single instanton results ( |10.12[ ) and 

dUD. 

Considering the quark-antiquark scattering we have to take into account that it is 
possible to treat an antiquark with velocity v 2 as a quark moving backward in time with 
velocity — v 2 . As a result, the scalar product of velocities changes the sign (vfvl) = —{v\v%) 
and the scattering angles are related as 



Then one gets 

d<Jgg{t) 

dt 



-.n r . 



Xqq -> «7T ~ Xqq 



coth 2 X F 2 (t) - l C -^F c (t)F s (t) + l ^-^F 2 (t) 



3 sinh x 



3 sinh x 



(10.28) 



. (10.29) 
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The second terms in ( 10.25 ) and ( 10.29| ) corresponding to the contribution of the C— odd 
amplitude has been missed in |53| because there only SU C (2) QCD is considered. 

The spin averaged total quark-quark cross section in the instanton-antiinstanton ap- 
proximation reads 



<J, 



gg 



dq 2 



FK q l) + \—F c { q l)F s {ql) 
6 s 



(10.30) 



which is constant in the high energy s limit. It is finite if the constraint instanton solution 
is used. In Eq. (10.30), the only term corresponding to the C = +1 exchange, Pomeron, 
survives at asymptotic, while the C = — 1 contribution, odderon, (second term in ( |10.30| )) 
is suppressed by the small factor ~ m 2 /s. The growing part of the total cross section arises 
as 

(10.31) 



(n c pi) 2 A(t)lns 



considering the inelastic quark - quark scattering in the instanton-antiinstanton background 

As it was discussed in detail in [16, 53 1, the model of the Pomeron based on nonpertur- 
bative gluon exchanges explains many properties of the diffractive scattering: the effective 
vector- like exchange ( |10.12 ), the additive quark rule and the main features of the total 
cross section ( |10.30| ), ( |10.31| ). 



It is important to note that the original Wilson loop ( 10.15 ) has essentially Minkowskian 
light-cone geometry whereas the instanton calculations of Wilson loop are performed in the 
Euclidean QCD. In analogy with the quark form factor the analytical continuation from 
Minkowski space to Euclidean one and v. v. becomes possible since the dependence of 
the Wilson loop on the total energy s and transverse momentum squared t is factorized in 
JlClrl ) and ( |j~Cl2l ). At high energy, the amplitude is s— independent both in the pertur- 
bative and nonperturbative cases. At the same time, the t— dependence of the amplitude is 
naturally expressed through the nonperturbative instanton field. Notice that in the origi- 
nal expression for the Wilson loop ( 10.15| ) , the nonlocal instanton correlator was integrated 
over both space-like and time-like separations x 2 corresponding to the distance between 
different points on the contour C qq , whereas the final expressions (|10.13|) , ( 110.271 ) depend 
on the space-like variable t. Thus we can proceed the formal calculations in Minkowski 
space and then make the Wick rotation k\ — > — k\, or, that is more natural from the point 
of view the instanton model, to perform formal manipulations in the Euclidean space and 
than make analytical continuation. 



11. Conclusions 

Besides the considerable progress in investigation of the role of nonperturbative QCD 
vacuum structure (in particular, of the instanton phenomena) in low and moderate energy 
domains of hadronic physics, nowadays there is a lack of understanding of their role in high 
energy processes which are intensively studied in modern experiments in particle physics. In 
this work we presented the results of analysis of the structure of nonperturbative corrections 
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in such important quantities as the quark form factor and the cross section of the diffractive 
quark-quark scattering at high energy. The quark scattering process was considered in 
background of QCD vacuum which is described within the instanton liquid model. The 
instanton contribution to the electromagnetic quark form factor is calculated in the large 
momentum transfer regime. We estimated analytically the weak field approximation for the 
instanton contribution and find the all-order result in the asymptotic regime. The latter was 
performed considering the Gaussian simulation of the constraint instanton profile function. 
We found that the leading contribution to high energy asymptotic behavior is provided 
by the lowest order terms. Although the latter result could be treated as an argument in 
favor of validity of the weak field approximation, the better work has to be done in this 
direction since the results for Gaussian profile and the instanton in the singular gauge may 
be different in general. 

The instanton effects are more important for phenomenology of the hadronic processes 
possessing two different energy scales. (For more detailed discussions see the works [pl||). 
One of such situations — quark-quark diffractive scattering — was considered in the last Sec- 
tion of the present work. This approach supports the model of the Pomeron as exchange by 
nonperturbative gluons interpolated by the classical instanton field. We have shown that 
the C— odd (odderon) amplitude is suppressed as 1/s compared to the C— even (Pomeron) 
one. In the case of diffractive scattering, the total center-of-mass energy s (hard charac- 
teristic scale) is large while the squared momentum transfer —t is small compared to the 
latter: —t <C s, but nevertheless larger than any IR scale. Besides this, the other cases of 
interest where the nonperturbative (including instanton induced) effects may be significant 
are the saturation in deep-inelastic scattering at small-x f73fl , and the transverse momen- 
tum distribution of vector bosons in the Drell-Yan process |4^] . The latter is one of the 
most important objects of the experimental investigations (in particular, in the context 
of searches for New Physics and Higgs bosons — at future LHC and Tevatron experiments 
|74|]), as well as theoretical studies of the predictive power of pQCD at various energy scales 
and the role of nonperturbative physics (see, e.g., (7!|] and references therein). 
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13. Appendix 

The evaluation of the path integrals in the weak-field (one-loop) approximation was per- 
formed by using the following expressions: The partial derivatives in n-dimensional space- 
time are presented as 

8 2 z = 2n8 z 2 + Az 2 2 z2 , d»d v = 2g^d z2 + A Zfl z u d 2 z2 . (13.1) 
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The scalar products of the scattering vectors ( |9.6[) in terms of the scattering angle x read 

v ji v l z t^ z f = z 2 coshx + ot sinh 2 x , 

z 2 = ( Vl a + v 2 T) 2 = a 2 + t 2 + 2<7Tcoshx. (13.2) 
To extract the universal cusp factor, we use the integrals: 

roc roo ^ 

da dr exp \— a(a +r +2o"Tcoshx)l = : — (13.3) 

Jo Jo 2a sum x 

/ do I dr ot exp \—a(a 2 + r 2 + 2ar cosh v)l = X c °t X ^ (13.4) 

Jo Jo 1 J 4a 2 sinh 2 x 

Therefore, for arbitrary functions Aj(z 2 ) we get: 



and 



oc 



and gets 



Af } (z 2 ) = (-) fc / daa fc e- Q * a A i (a) (13.5) 



oo />oo 
oo roo 

da / dr z 2 A"(z 2 ) = — ^— A(0) , (13.6) 
o Jo 2sinh X KJi V ; 

^daTdr arA"(z 2 ) = - 1 - XC f X A(0). 
o Jo 4 sinh x 

In the all-order instanton calculations, the Gaussian integrals over zq are taken by 

'dV(4) n e-i =n! _J_, (13.7) 

where 

pn = n! 
s s!(ra-2s)!2 2 * 



and we use the binomial formula 

n 

zi + zl) n = Y^C n k zfzT- k) . (13. 
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Figure 1: The notations of the quark momenta for the (a) quark form factor, and (b) quark-quark 
scattering. 
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(a) (b) (c) 




Figure 2: The total cusp-dependent part of the Wilson loop integral for the quark form factor (a); 
the leading order contributions of the perturbative (b) and nonperturbative (single-instanton) (c) 
fields; (d) the all-order single instanton result; (c) the exponentiation of the single instanton result. 
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Figure 3: The asymptotic behaviour of the quark form factor is shown as the function of the di- 
mensionless variable q 2 — Q 2 /A 2 , up to terms O (lnlnq 2 ). The long-dash presents the contribution 
of one loop perturbative terms; the solid line represents the total form factor including the instanton 
induced part, Eq. (8.21). For comparison, the leading (~ In q 2 In In q 2 ) perturbative contribution is 
shown separately — the short-dash line. 
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Figure 5: Lowest orders instanton contributions to the Wilson integral with spin factors. The 
leading term (a\a\) (7) is the dashed line, next-to-leading terms (ajo^) (7) and (cuf of) (7) are the 
dotted and dash-dotted line, correspondingly. The sum of these contributions is the solid line. 
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Figure 6: Lowest orders instanton contributions to the Wilson integral without spin factors. The 
leading term (a\a\) (7) is dashed line, next to leading terms (aja^) (7) and (afaf) (7) are dotted 
and dash-dotted line, correspondingly. The sum of these contributions is the solid line. 
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(a) (b) (c) 

Figure 7: Schematic representation of partial summation of the instanton contributions. 
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